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ABSTRACT 

The cosmic history of the growth of supermassive black holes in galactic centers parallels 
that of star-formation in the Universe. However, an important fraction of this growth occurs 
inconspicuously in obscured objects, where ultraviolet/optical/near-infrared emission is heav¬ 
ily obscured by dust. Since the X-ray flux is less attenuated, a high X-ray-to-optical flux 
ratio ( fxl fo ) is expected to be an efficient tool to find out these obscured accreting sources. 
We explore here via optical spectroscopy. X-ray spectroscopy and infrared photometry the 
most extreme cases of this population (those with f x jf 0 > 50, EXO50 sources hereafter), 
using a well defined sample of 7 X-ray sources extracted from the 2XMM catalogue. Five 
EXO50 sources (~ 70 percent of the sample) in the bright flux regime explored by our sur- 
vey (/(2-iofev) > 1.5 x 10 -13 erg cm -2 s -1 ) are associated with obscured AGN (Nh > 10 22 
cm -2 ), spanning a redshift range between 0.75 and 1 and characterised by 2-10 keV intrinsic 
luminosities in the QSO regime (e.g. well in excess to 10 44 erg s -1 ). We did not find com¬ 
pelling evidence of Compton Thick AGN. Overall the EXO50 Type 2 QSOs do not seem to 
be different from standard X-ray selected Type 2 QSOs in terms of nuclear absorption; a very 
high AGN/host galaxy ratio seems to play a major role in explaining their extreme proper¬ 
ties. Interestingly 3 out of 5 EXO50 Type 2 QSO objects can be classified as Extreme Dust 
Obscured Galaxies (EDOGs, fiApmlfR > 2000), suggesting that a very high AGN/host ratios 
(along with the large amount of dust absorption) could be the natural explanation also for 
a part of the EDOG population. The remaining 2 EXO50 sources are classified as BL Lac 
objects, having rather extreme properties, and which are good candidates for TeV emission. 

Key words: galaxies; active - galaxies: neclei - BL Lacertae objects; general - Quasars; 
general - X-rays; galaxies 


1 INTRODUCTION 

The study of high-z obscured quasars (Type 2 QSOs: the high 
luminosity counterpart of Seyfert 2 galaxies) is one of the hot 
topics of current extragalactic astronomy. Their observed distri¬ 
butions (e.g. redshift, flux and absorption) and physical proper¬ 
ties (e.g. bolometric luminosity, black hole mass, accretion rate), 
compared to those of unobscured QSOs, are key parameters to 
unders tand and t o test the orientation-based Unified Schemes 
d AntonuccTl 1 1 9931) for Active Galacic Nuclei (AGN) and to con- 
stra in the contribution of QSOs to the X- ray Background (XRB: 
e.g. iGilli et al] 120071 : iTreister et alj|2009l) . Even more important, 
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there are now increasing evidence that the formation and growth 
of galaxies and their nuclear supermassive black holes (SMBHs, 
M BH > 10 6 M q ) are intimately related; feedback from actively 
accreting SMBH, AGN, must play a fundamental role in regu¬ 
lating bot h star formation and accretion throughout galaxy’s his- 


2005 

, Croton et al. 200q^ Hopkins et al. 

20081 Menci et al. 12008S, 

King 

2010|, Faucher-Giguere & Quataert 

2013). Although the na- 


ture of this relationship is still poorly understood, there are hints 
that most of the SMBH accretion takes place during an obscured 
quasar phase. The infrared (IR), optical and X-ray spectral analy¬ 
sis of a large sample of obscured quasars probably represents one 
of the best methods to test a possible c oevol u tion between mas- 
sive galaxies and AGN activity (see e.g. lHopkins et al.ll200d) and 
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to investigate if and how the AGN feedback can affect the galaxy 
evolution (see e.g. Bower et al] |2006l ). Indeed, in these sources the 
properties of the host galaxy can be studied using the optical and 
near-IR data, where the absorbed AGN is supposed to contribute 
marginally, while the AGN properties can be investigated using 
the X-ray and mid-IR data where, conversely, the AGN emission is 
dominant (see lBallo et al.lEoi4l and references therein). Obscured 
QSOs, a rare class of objects, are thus expected to have large val¬ 
ues of the X-ray-to-optical flux ratio. Furthermore, since the dust 
extinction increases in the ultraviolet (UV) while the X-ray absorp¬ 
tion goes in the opposite direction, i.e. strongly decreasing going 
towards the high ener gies, a redshift dependence is also expected 
(~ (1 + z) 3 ' 6 : lFiore et al.l2003l) . However, in spite of the outstanding 
pro gress, o btained by using medium and deep X-ray surveys (see 
e.g. Brusa et al ■201 Qand references therein), the weakness, both in 
the X-ray and in the optical band, of the selected sources with high 
values of the X-ray-to-optical flux ratio usually p revents from a de¬ 
tailed analysis of the ind ividual objects (see e.g. iKoekemoer et~akl 
l2004ICivano et al.ll2005l) . 

To explore the most extreme examples of obscured QSOs 
in the bright flux regime, we have started a project focused on 
the very “high” X-ray-to-optical flux ratio, f x lf 0 (see Eq. 1), 
population, namely those sources with fxlfo >50 (more than 
15 times the average values of unobscured broad line AGN, e.g . 
ICacci anig a et al. l2004 De lla Ceca et al. H2004l .[ci vano et al.l201 2l ) 

and with f x > 10~ lj erg cnT 2 s~' . In this way we should efficiently 
select the best candidates to be obscured QSOs; at the same time, 
the source brightness ensures that the optical spectroscopic identi¬ 
fication can be achieved Q for sources with reasonably good qual¬ 
ity X-ray data (>few hundred counts) to carry out a reliable X-ray 
spectral analysis. Hereafter we will call these sources EXO50 - for 
“extreme X-ray to optical flux ratio”. One example of such extreme 
obscured objects is XBS J021642.3-043553 with f x /f 0 ~ 200, for 
which the presence of a Type 2 Q SO at z ~ 2 was spec troscop- 
ically confirmed with VLT/FORS dSevergnini et aljl200 6). Using 
the redshift information and the spectral energy distribution (SED), 
we estimated a stellar mass of ~ 10 11 M G for the host galaxy, which 
supports the strong link between high-redshift massive galaxies and 
powerful obscured high-redshift QSOs. In spite of the importance 
of obscured QSOs in the cosmological context, only less than a 
dozen of such extreme f x /f 0 sources have been found and stud- 


ied so far with good X 

ray and optical data (see e.g. Gandhi et al. 

2006 

. Severgnini et al. 

2006, Del Moro et al. 2009L ICarmfisi et al. 

2009 

. Brusa et al.l 2010 

Brusa et al. 20151, Perna et al. 2014). 


Here we discuss the first results obtained from this project on 
a small, but statistically complete and representative, sub-sample of 
7 EXO50 objects with f x > 1.5 x 10 -13 erg cnT 2 s -1 . It is worth 
mentioning that other two interesting (and rare) classes of extra- 
galactic sources are expected to show up in the high f x lf 0 domain 
explored here, namely BL Lac objects and high redshift (i.e. z> 0.6 
in order to have the 4000 A break shifted at longer wavelengths 
than the R band filter) clusters of galaxies. Clusters of galaxies, 
however, are not expected to be represented in our sample since we 
have restricted our search only to the point-like X-ray sources, thus 
strongly minimising their possible selection. 

This paper is organized as follows: in Section 2 we discuss 
the strategy used to define a statistically representative sample of 
sources with f x !f 0 >50. In Section 3 we present the data accumu¬ 


lated so far (in the radio, IR, optical, and X-ray energy ranges) on 
the 7 EXO50 objects discussed here, while their properties, source 
by source, are discussed in Section 4. In Section 5 we compare the 
broad band, from IR to X-rays, properties of these EXO50 sources 
with other similar objects found in the literature. Summary and 
conclusions are presented in Section 6. Throughout this paper, we 
consider the cosmological model with (H„, Qm,Qi)=( 70,0.3,0.7); 
results from other papers have been rescaled to this cosmological 
framework. All the optical/IR magnitudes reported here are in the 
Vega system. Unless differently quoted, X-ray luminosities are in¬ 
trinsic (i.e. unabsorbed) luminosities in the rest-frame 2-10 keV 
energy range. In this paper we use the term Type 1 and Type 2 
AGN as broad or narrow line (FWHM of the permitted emission 
lines < 1500 km s _1 ) AGN, irrespective to their intrinsic X-ray 
luminosity, while we use the term Type 1 QSO and Type 2 QSO 
(or obscured QSO) as broad or narrow AGN with an intrinsic 2- 
10 keV luminosity in excess to 10 44 erg s _1 . Finally, errors are at 
90% confidence level for the X-ray spectral parameters derived us¬ 
ing XSPEC (as usually done in X-ray astronomy) and 68% for all 
the other quantities. 


2 THE EXO50 SAMPLE 


In obscured QSOs the nuclear UV/optical emission is suppressed 
by dust obscuration (thus leaving only the galaxy component vis¬ 
ible), whilst the nuclear X-ray flux (e.g. in the 2-10 keV energy 
range), even if exhibiting signat ures of pho toelectric absorption, 
is less attenuated (see e.g. iFiore et aPl2003l) . Obscured QSOs are 
therefore expected to display large values of X-ray to optical flux 
ratio, defined here as: 

f x lfo = A/(A/1 x 2.15 x 10~ 9 x 10- 0 ' 4 **) (1) 


where f x refers to the observed 2-10 keV flux (corrected for Galac¬ 
tic absorption), R is the observed opt ical magnitude in the R band 
(T = 6410 A ) and A A - 1568 A (see lFukugita et all 19951) . 

We define EXO50 the sources with /*/ f B > 50; these are a 
rare class of X-r ay emitting objects, as ty p ically AGN have f r / f„ 

1 1 1 i (see e.g. Caccianiga et al. I 120041 . iDella Ceca et al. 11 20041 . 
ICivano e t all201 2h. For instance, in the XMM-Newton Br ight Sur¬ 
vey (XBS. lDella Ceca et al. 1120041 [Caccianiga et alj|2008h . a com- 
plete sample of bright (/( 0 . 5 - 4 . 5 keV) ^ 7x 10~ 14 erg cnT 2 s _1 ) X-ray 
selected sources almost completely identified (spectroscopic iden¬ 
tification level ~ 98%), EXO50 objects represents only ~ 0.5% of 
the source population. 

To construct our sample of EXO50 sources in the bright flux 
regime we used one of the largest well defin ed and com p lete X- 
ray source sample derived so far (disc ussed inlMate os et alJ l2008l) . 
based on the 2XMM source catalog 0 l lWatson et all200& . 

First we have considered the following selection criteria: 

a) “point-like” X-ray sources (parameters EPJ3XTENT and 
EP_EXTENT_ML in the source catalogue equal to 0), in order to 
minimise the presence of clusters of galaxies which is another pos¬ 
sible class of high f x l f 0 sources; 

b) north of -10 deg declination, in order to be accessible to 
optical investigation from the Italian Telescopio Nazionale Galileo 
(TNG), from the Large Binocular Telescope (LBT) and from the 
Spanish Gran Telescopio Canarias (GTC); 


1 At the chosen X-ray flux limit of f x ~ 10 13 erg cm 2 s 1 a source with 

fx!fo ~ 50 (300) has an optical R magnitude of ~23 (~25). 2 See http://xmmssc-www.star.le.ac.uk/Catalogue/xcat_public_2XMM.html 
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c) sources selected in the 2-10 keV energy band with a/ 2 _io kev 

> 10~ 13 erg cur 2 s _1 ; 

d) only serendipitous sources (e.g. sources that are not related 
to nearby galaxies or to the XMM-Newton pointing) have been 
taken into consideration. We have also excluded the serendipitous 
sources already classified as non AGN. 

_Starting from a source list of 9431 sources (see lMateos et al.l 

l2008h . these first selection criteria provide us with a list of about 
600 X-ray emitting objects. 

The second step was to select the X-ray sources with 
fjfo >50, using: 

i) optical source archives (e.g. APM, Sloan Digital Sky Survey 
- SDSS, etc); 

ii) optical imaging data from archives (STScI Digitized Sky 
Survey, ESO archive, etc..); 

iii) imaging data from our dedicated observing programs (e.g. 
TNG imaging). 

We have looked at and visually inspected all the selected X- 
ray sources for possi bl e optical co unterparts, using a search radius 
of 4 arcsec (see lCaccianiga et~ai]|2008l) . and derived their R mag¬ 
nitude. For all the sources reported and discussed here (see Table 
1), the offset between the X-ray and the optical position, derived 
a posteriori , is below 2.1 arcsec, fu lly consistent with the results 
obtained bv lCaccianiga et all |2008| ) from the analysis of the X-ray 
to optical offset for the X-ray sources belonging to the XBS. The 
accumulated magnitudes, combined with the observed 2-10 keV 
fluxes, have been used to compute the f x lf 0 ratio and thus to select 
EXO50 sources. 

In this paper we consider a first complete sample of EXO50 
sources having f x > 1.5 x 10 -13 erg cm -2 s -1 ; this sample is com¬ 
posed by the 7 EXO50 objects reported in Table 1. In Figure 1 
we show the f x lfo distribution of the sample composed by all 
the serendipitous point-like sources in the covered area having f x 

> 1.5 x 10 -13 erg cm~ 2 s _1 (i.e. the starting sample, 310 objects). In 
the same figure we show how extreme the 7 sources are with respect 
to the starting sample; the EXO50 sources discussed here represent 
only ~ 2% of the population considered. We stress that according 
to the selection criteria discussed above and the optical material 
at our disposal these 7 objects are the only EXO50 sources, north 
of -10 deg declination, wit h f x > 1.5 x 1 CL 13 erg cnr 2 s _1 in the 
sky survey area defined in lMateos et al.l ( 12008 3; the sky coverage 
investigated to find out these 7 objects is about 60.4 sq deg. 

As detailed in the following sections, out of the 7 EXO50 
sources 2 are BL Lac objects, 3 are confirmed Type 2 QSO and 
2 remain spectroscopically unidentified, although their broad band 
properties strongly suggest an ob scured Q SO nature . Us ing the 
JAVA software GALAXYCOUNT ( lEllis & Bland-Hawthoml2007l) 
and considering the optical magnitudes of the proposed optical 
counterparts (see Table 1) we estimate less than 0.3 normal galaxies 
by chance in the total area covered from the 7 error boxes, assuring 
us about the reliability of the spectroscopic identification proposed 
here. 


3 OBSERVATIONS AND DATA REDUCTION 
3.1 Optical imaging 

In Figure 2 (left panels) we present the optical finding charts of the 
7 sources discussed here (lxl arcmin) produced using data from 
dedicated TNG observing runs (R filter) or from the SDSS (r filter). 
TNG raw images were reduced following the standard CCD 
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Figure 1 . Histogram of the f x /f 0 distribution of the starting sample with f x 
> 1.5 X 10 13 erg cm~ 2 s _1 (310 objects). We also show as filled histogram 
the f x /f 0 distribution of the 7 EXO50 sources discussed in this paper. 


reduction process using IRAF2A4. All images were de-biased and 
corrected for flat-field effects. We combined the different exposures 
for each source and the resulting images were flux calibrated using 
standard stars and the standard extinction curve for the observa¬ 
tory. Finally, we did an astrometric calibration of each image; the 
astrometric accuracy of our imaging is typically better than ~ 0.3 
arcsec. The R magnitudes were estimated with IRAF routines fol¬ 
lowing the standard procedure of aperture photometry centered at 
the source position. 


3.2 Optical spectroscopy and identification 

The optical spectra for 6 of our sources were taken at GTC dur¬ 
ing 3 different observing runs (see Table 2) using the R500R 
grism on the OSIRIS instrument 0. For the reduction we used 
the standard IRAF long-slit package following the standard steps. 
The resulting spectra (see Figure 2 middle panels, in arbitrary 
units) give us spectroscopic information (classification and red- 
shift) for 4 out of the 6 GTC sources. For the remaining 2 sources 
(2XMMJ 100038.9+050955 and 2XMMJ135055.7+642857) the 
spectra are inconclusive and do not allow us to derive any rele¬ 
vant information. The source 2XMMJ121026.5+392908 is a well 
known object and its redshift and classification have been taken 
from the literature. In total we have spectroscopically identified 
(using the GTC data as well as data from the literature) 5 EXO50 
objects; two sources remain still spectroscopically unidentified, al¬ 
though the data at other wavelengths discussed in section 4 and 5 
strongly suggest an obscured QSO nature. 


3 See http://www.gtc.iac.es/instruments/osiris/osiris.php 
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Figure 2. Left panels : Optical finding chart (lxl arcmin) for the 7 EXO50 objects reported in this paper produced using material from dedicated TNG 
observing runs (R filter; sources: 2XMM022256.9-024258, 2XMM 100038.9+050955, 2XMM 123204.9+215254, 2XMM 135055.7+642857) or from SDSS 
(r filter; sources: 2XMM121026.5+392908, 2XMM121134.2+390054, 2XMM143623.8+631726). The circle is 4 arcsec radius around the X-ray position. 
Middle panels : The optical spectra of the EXO50 objects observed at the GTC. The solid black line in the spectroscopically identified objects represents the 
underline galaxy component used to reproduce the optical spectrum, except for 2XMMJ121134.2+390054, for which we show the galaxy component (blue 
line), the AGN component (continuum plus broad emission lines, red line) and the combined component (galaxy plus AGN; black line). The spectral lines used 
to classify the objects are also marked. 2XMMJ121026.5+392908 is a well known object and its redshift and classification has been taken from the literature. 
Right panels: X-ray data (observer frame) and residuals for the EXO50 objects discussed here. All the objects are well fitted with an absorbed power-law 
model. Red filled triangles: pn data. Black open circles: MOS data. 
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Table 1. The complete sample of EXO50 objects discussed in this paper 


n 

2XMMName 

fx 

R 

fx/fo 

ID 

z 

N a „ 

r fl 

L x 

fr 



10 13 cgs 





10 22 cm - 2 


10 45 cgs 

mJy 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(id 

i 

J022256.9-024258 

3.50 

22.1+0.4 

72 

AGN2 

1.004 

7 5 +3-2 
' -2.9 

00 

o 

1 + 

© p 

1.9 

429+14 

2 

J 100038.9+050955 

1.63 

23.6+0.4 

138 

- 

- 

>0.4* 

1.9 C 

- 

- 

3 

J121026.5+392908 

36.1 

19.16 rf +0.02 

50 

BLLac 

0.617 

0.06 + °“; 

2.23 + °'°! 

6.4 

19.0+0.7 

4 

J121134.2+390054 2 

8.2 

20.77 rf +0.05 

50 

BLLac 

0.890 

0.19+°® 

2 - 21+ ool 

3.7 

10.6+0.6 

5 

J123204.9+215254 e 

10.0 

23.9+0.4 

1118 

AGN2 

0.763 

3.61 -’ll 

1.3 Up!:? 

1.9 


6 

J135055.7+642857 

1.5 

25.0+0.8 

458 

- 

- 

>0.1* 

2 -° 2 io 


183.5+5.5 

7 

J143623.8+631726 

2.53 

22.2 rf ±0.2 

55 

AGN2 

0.893 

1 46 +0 - 39 
-0.24 

0.84 

- 


Columns are as follows: 1) Number used to mark the object in the plots shown in Section 5; 2) source name; 3) X-ray flux in the 2-10 keV energy band (MOS 
normalisation) corrected for Galactic absorption in units of 10“ 13 erg cm -2 s _1 ; 4) R band magnitude; 5) X-ray to optical flux ratio; 6) optical spectroscopic 
classification; 7) redshift; 8) intrinsic absorption and 90% confidence interval in units of 10 22 cm -2 from fits to X-ray spectra; 9) power law photon index, T, 
in the X-ray energy range, and 90% confidence interval; 10) intrinsic X-ray luminosity in the 2-10 keV rest frame energy range, in units of 10 45 erg s -1 ; 11) 
radio flux at 1.4 GHz in mJy. Notes: a - For all the sources the best fit spectral model in the X-ray energy band is a simple absorbed power-law model; b - The 
90% lower limit onN// obtained assuming z=0; c - This parameter has been fixed; d - For these sources the R band magnitude has been obtained from the 
SDSS r ’’model” magnitudes assuming r-R = 0.27 ( tsdss ab = Rvega + rAfl(Vega) with r / 4fi(Vega)=0.27); e - For these two sources we propose here a 
different redshift with respect to that already reported in the literature. See Section 4 for details. 


Table 2. Optical follow-up spectroscopic observations 


2XMMName 

(1) 

Programme 

(2) 

Instrument/Grism 

(3) 

Exposure (s) 
(4) 

Slit Width (") 
(5) 

Dates 

(6) 

J022256.9-024258 

GTC18-10B 

OSIRIS/R500R 

6x1100 

1.2 

2-Sep-2010 (4), 19-Sep-2010 (2) 

J100038.9+050955 

GTC44-11A 

OSIRIS/R500R 

4x1800 

1.0 

10-Apr-2011 

J121134.2+390054 

GTC18-10B 

OSIRIS/R500R 

2x600 

1.2 

27-Jan-2011 

J123204.9+215254 

GTC44-11A 

OSIRIS/R500R 

2x1800 

1.0 

10-Apr-2011 (1), 11-Apr-2011 (1) 

J135055.7+642857 

GTC44-11A 

OSIRIS/R500R 

3x1800 

1.0 

10-Apr-2011 

J143623.8+631726 

GTCMULTIPLE2-09A 

OSIRIS/R500R 

7x1500 

1.2 

17-Apr-2009 (1), 19-Apr-2009 (3), 25-May-2009 (3) 


Columns are as follows: 1) Source name; 2) program ID; 3) instrument/grism used 4) number of exposures and exposure times in seconds; 5) slit width in 
arcsec; 6) dates of the observations. In parenthesis we have indicated the number of exposures. 


3.3 WISE All-Sky Survey data 

The Wide-field Infrared Survey Explorer (WISE; IWrightetaD 
l2dIob has recently carried out an all sky survey in the medium in¬ 
frared regime, detecting hundreds of millions of objects. The es¬ 
timated 5cr point sources sensitivities (in unconfused regions) in 
the four observed channels (Wl=3.4 /im, W2=4.6 /rm, W3=12 pm 
and W4=22 pm) are better than 0.08, 0.11, 1 and 6 mJy, while 
the angular resolution (FWHM) are 6.1, 6.4, 6.5, 12 arcsec re¬ 
spectively; the astrometric precision for high signal-to-noise ratio 
(SNR) sources is better than 0.15 arcsec. To add the WISE infor¬ 
mation to our study we have cross-correlated our EXO50 sample 
with the WISE All-Sky source catalog [] using a positional toler¬ 
ance, from the X-ray source position, equal to 4 arcsec . All our 
sources have been detected by WISE (see below) and a posteri¬ 
ori the offset between the WISE and the optical position of our 
sources is below 2 arcsec. The results of this cross-correlation are 
reported in Table 3. We found a single WISE counterpart for all 
the 7 sources discussed here. The detections reported in Table 3 
have a significance greater than la in all bands, with the excep¬ 
tion of 2XMMJ100038.9+050955, 2XMMJ123204.9+215254 and 
2XMMJ143623.8+631726, detected in the W4 band at 6.1cr, 3.6cr 


4 We use here the public available All-Sky Data Release 
that covers >99% of the sky (March 2012 release; see 
http://wise2.ipac.caltech.edu/docs/release/allsky /1 


and 3.9cr respectively. Two sources, 2XMMJ 121026.5+392908 
and 2XMMJ121134.2+390054, are not detected in the W3 and W4 
bands. 

Using a simulated sample with random positions of ~ 1600 
object we estimated that the probability to have a random WISE 
source inside a circle of 2 arcsec radius is ~ 1% (i.e. 0.07 WISE 
sources expected by chance in the 7 error circles investigated here), 
implying that all the detected WISE counterparts are very likely 
associated with the high f x /f a sources. 

Flux densities at 3.4 pm, 4.6 pm, 12 gm and 22 pm have been 
computed from the magnitudes reported in Table 3 by assuming 
the magnitude zero points of the Vega syste m correspo n ding to a 
power-law spectrum (/„ oc y“) with a = 1 (see lWright et al.l2010h . 
The differences in the computed flux densities expected using flux 
correction factors that correspond to a A a = ±1 are lower than 0.8 
per cent, 0.6 per cent, 6 per ce nt, and 0. 7 per cent in the Wl, W2, 
W3 and W4 band, respectively l l Wright e~ all20ld) . 

3.4 XMM-Newton spectroscopy 

In Table 4 we report details for the XMM-Newton data used for the 
X-ray spectral analysis of each source discussed here. The XMM- 
Newton data were cleaned and processed with the XMM-Newton 
Science Analysis Software (SAS) and were analyzed w ith standard 
software packages (FTOOLS; XSPEC. lAmaud IT996I ). Event files 
produced from the pipeline were filtered from high-background 
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Table 3. Infrared (WISE) data for the EXO50 sources discussed here. 


n 

(1) 

2XMMName 

(2) 

fjfo 

(3) 

ID(z) 

(4) 

W1 

(5) 

W2 

(6) 

W3 

(7) 

W4 

(8) 

Log (vL v ) 
(12.3/tm) 
(9) 

i 

J022256.9-024258 

72 

AGN2(z= 1.004) 

16.02 ± 0.06 

14.89 ± 0.07 

11.08 + 0.09 

8.07 + 0.15 

45.61 

2 

J100038.9+050955 

138 

- 

16.22 + 0.08 

14.29 ± 0.05 

10.33 + 0.06 

7.97 + 0.18 

- 

3 

J121026.5+392908 

50 

BLLAC(z=0.617) 

14.86 ± 0.04 

14.59 ± 0.06 

> 12.34 

> 8.66 

- 

4 

J121134.2+390054 

50 

BLLAC(z=0.890) 

15.21+0.04 

15.03 ± 0.08 

> 12.50 

> 9.45 

- 

5 

J123204.9+215254 

1118 

AGN2(z=0.763) 

15.35 + 0.05 

13.84 ± 0.04 

10.19 + 0.05 

8.48 + 0.30 

45.10 

6 

J135055.7+642857 

458 

- 

14.90 ± 0.03 

13.85 + 0.03 

10.46 ± 0.04 

8.47 + 0.15 

- 

7 

J143623.8+631726 

55 

AGN2(z=0.893) 

15.09 + 0.03 

14.05 ± 0.04 

11.36 + 0.09 

8.91+0.28 

45.10 


Columns are as follows: 1) Number used to mark the object in the plots shown in Section 5; 2) source name; 3) X-ray to optical flux ratio; 4) optical 
spectroscopic classification and redshift; 5) W1 magnitude (3.4pm) and I <r error; 6) W2 magnitude (4.6 pm) and lcr error; 7) W3 magnitude (12 pm) and lcr 
error. The magnitude lower limits correspond to the 95% confidence level as reported in the WISE catalog; 8) W4 magnitude (22 pniI and lcr error. The 
magnitude lower limits correspond to the 95% confidence level as reported in the WISE catalog; 9) Log of the rest frame 12.3 /cm luminosities 

(vi 2 . 3 um x Li 2 . 3 /j m in erg s -1 ) computed as described in Section 3.3. 


Table 4. EPIC XMM-Newton observation details. 


Name OBSID 

(1) (2) 

Nh Instr. 

(3) (4) 

Filter 

(5) 

Net Exp. Time 
(6) 

Net Counts 
(7) 

J022256.9-024258 0037981601 

2.34 MOS1&2 

thin 

26.1 

265 


pn 

thin 

7.9 

127 

J100038.9+050955 0204791101 

2.41 MOS1 

med 

16.1 

27 


pn 

med 

12.5 

63 

J121026.5+392908 0112830501, 2.00 MOS1&2 

med 

161.5 

119876 

0112830201 

pn 

med 

68.9 

175724 

J121134.2+390054 0112190201 

1.84 MOS1&2 

med 

27.8 

2311 


pn 

med 

8.9 

2488 


J 123204.9+215254 0112650301 1.80 MOS1 

thin 

19.3 

248 

MOS2 

med 

19.5 

250 

pn 

thin 

13.1 

447 


J135055.7+642857 0147540101 

1.70 MOS1&2 

med 

44.6 

873 


pn 

med 

18.9 

1064 

J143623.8+631726 0204400301 

1.37 MOS1&2 

med 

46.7 

862 


pn 

med 

12.7 

648 


Col. (1): 2XMM source name; Col. (2): OBSID of the XMM-Newton observation; Col. (3): Galactic absorbing column density along the line of sight in units 
of 10 20 cnT 2 ; Col. (4): EPIC instrument; Col. (5): EPIC filter; Cols (6): Exposure time after removing high-background intervals, in units of ksec; Col. (7): 

Net counts in the energy range 0.3 - lOkeV. 


time intervals and only events corresponding to patterns 0-12 for 
MOS and 0-4 for pn were used. All spectra were accumulated from 
a circular extraction region with a radius of 20-30 arcsec, depend¬ 
ing on the source off-axis distance. Background counts were accu¬ 
mulated in nearby circular source free regions, using an area usu¬ 
ally about a factor of 4 larger than the one used to extract the source 
counts. The X-ray spectra usually cover the 0.3-10 keV energy 
range; the total (MOSl+MOS2+pn) net (background subtracted) 
counts range from ~ 10 2 to ~ 3 x 10 5 counts. The ancillary re¬ 
sponse matrix and the detector response matrix were created by the 
XMM-SAS tasks arfgen and rmfgen at each source position in the 
EPIC detectors. To improve the statistics, the MOS1 and MOS2 
spectra obtained by using the same filter were combined a poste¬ 
riori by using the FTOOLS task mathpha; in this case ancillary 
and detector response matrices for the MOS1 and MOS2 detectors 


for each source were combined by using the tasks addrmf and ad- 
darf. For all the sources but 2XMMJ100038.9+050955 we grouped 
the spectra in bins containing more than 20 (source+background) 
counts and used the % 2 minimization technique; in the case of 
2XMMJ100038.9+050955 we grouped the spectra in bins contain¬ 
ing 10 (source+background) counts and use the Cash statistics. 
In the case of 2XMMJ121026.5+392908 two XMM observations 
were used; pn (MOS) data from the two independent data set were 
combined together and the ancillary and detector response matrices 
were created using the same procedures quoted above. We fitted 
pn and MOS spectra simultaneously in the 0.3-10 keV band, tying 
together all pn and MOS parameters except for a relative normal¬ 
ization, which accounts fo r the differen ces between pn and MOS 
flux calibrations (see lMateos et al.ll2009l) . In the following, derived 
fluxes and luminosities refer to the MOS instrument. 
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For the spectral modelling we considered a simple ab¬ 
sorbed power-law model that takes into account both the Galac- 
tic hydrogen column de nsity along the line of sight (from 
iDickev & Lockmanll990l) and a possible intrinsic absorption at the 
so urce redshift (abundances relative to the Solar one as reported 
in IWilms et alj|200ch . In the X-ray spectral modelling we made 
use, when available, of the redshifts obtained from the optical spec¬ 
troscopy. All the X-ray spectra are well fitted by this simple model; 
the results are reported in Table 1 along with the corresponding flux 
(corrected for Galactic absorption) and intrinsic luminosity (i.e. 
corrected for both Galactic and intrinsic absorption) in the stan¬ 
dard 2-10 keV energy band. The X-ray spectra are shown in Figure 
2 (right panels). In the case of 2XMMJ 100038.9+050955, the spec¬ 
tral quality does not allow us to constrain the power-law photon in¬ 
dex and the intrinsic absorption at the same time, so we have fixed 
the power-la w photon index T to 1.9, a common value for unab- 
sorbed AGN i lMateos et a] .1120fol . ICorral et all 201 il iLanzuisi et al.l 
120131) . 


4 RESULTS 

Before discussing the properties of each single object in this sam¬ 
ple we summarise here their main properties. From optical spec¬ 
troscopy, out of the 7 EXO50 sources 2 are classified as BL Lac 
objects, 3 are classified as Type 2 QSO and 2 remain unidentified. 
The 3 sources classified as Type 2 QSO are in the redshift range 
0.7-1 and are characterised by an intrinsic X-ray absorbing col¬ 
umn density, N H , in the range between 1.5 and 8 xlO 22 cnr 2 . For 
the two unidentified objects the X-ray analysis provides lower lim¬ 
its to the intrinsic Nh of 4 x 10 21 cnT 2 and 10 21 cnT 2 , respectively. 
There are no Compton Thick AGN ( N H in excess to ~ 10 24 cnr 2 ) 
amongst the 3 Type 2 QSO, neither is suspected their presence in 
the 2 still spectroscopically unidentified sources, since their X-ray 
spectra are at odds with that usually observed in Compton Thick 
AGN (e.g. the presence of a prominent Iron line at 6.4 keV, rest 
frame, or a very flat X-ray spectrum). The 2 BL Lac objects are at 
z=0.62 and z=0.89 and, as it will be discussed in Section 5.7, they 
are rather extreme in their SED. 

For the three EXO50 sources spectroscopically identified as 
Type 2 QSOs (see below), we report in Table 3 the rest frame 
12.3 fjm luminosity obtained by interpolating the observed lumi¬ 
nosity (i.e. not corrected for reddening) in the W3 (observed frame, 
12 pm) and W4 (observed frame, 22 pm) WISE bands and us¬ 
ing the measured spectral index between 12 pm and 22 pm; the 
same spectra l index has been used to evaluate the K-correction. 
iGandhi et al.l J2009h . studying a sample of Type 1 and Type 2 
Seyferts, showed that the observed luminosity around 12.3 pm 
rest frame 612 . 3 ^ x should represent an accurate proxy 

for the AGN intrinsic power; in powerful AGN, the contribution 
expected from the host_galaxy at these wavelengths is marginal 
(< 10 %: lBallo et al.|[2oT4T) . Overall the derived infrared luminosi¬ 
ties are in the range typical of Luminous Infrared Galaxies (LIRG: 
L lR > 10 n L G ) and Ultra Luminous Infrared Galaxies (ULIRG: 
Lir > 10 12 L o , see e.g. lSanders & MirabeJl996h . Assuming abolo- 
metric correction factor to the 12.3 pm lumi nosity of ~ 10.8, appro¬ 
priate for high luminosity AGN (see iBallo et al.i[20l4h the implied 
bolometric luminosities for the three spectroscopically identified 
Type 2 QSO are in the range 1.5 - 4.5 x 10 46 erg s -1 . 

• 2XMMJ022256.9-024258; fjf 0 =72 

A single source (R=22.1) is evident in the optical finding chart 



Figure 3. Radio (FIRST) contour plots superimposed to a CCD R image 
(taken at the TNG) of 2XMMJ022256.9—024258; the red hatch indicates 
the object responsible for the X-ray emission. The radio contours start from 
0.2 mjy with a multiplicative step of ~ 1.58. The large blue bar at the bottom 
is 1' wide: at the redshift of the source T corresponds to ~ 0.48 Mpc. North 
is upwards. East is on the left. 


within the X-ray error circle (see Figure 2). The optical spec¬ 
trum of this object shows several narrow (observed FWHM < 
1000 - 1500 km s -1 ) emission lines that we associate with 
MgIL12798, [NeV]T3346,3426, [OII]T3728, [NeIII]T3869, H f , H*, 
Hy, [OIIIj/14364 and HeLl4471 at z = 1.004. The permitted emis¬ 
sion lines (Hy and MgII/12798) have widths similar to the forbidden 
ones. The Hy has a FWHM of 550-600 km/s that is close to the in¬ 
strumental resolution. The MgII/12798 has a strongly asymmetric 
profile, which is difficult to analyse with the present data, being the 
line at the border of the sampled wavelength range; the FWHM of 
this line is likely below 1500 km/s. It is worth noting the presence 
in this object of the [NeV]T3426 line, a reliable signature of nu¬ 
clear activity sin ce it can not be produced in starburst /starforming 
galaxies (see e.g. iGilli et alJl20HllMignoli et alj|2013l) . The X-ray 
spectrum is well described by an absorbed power-law model with 
an Nu of ~ 7.5 x 10 22 cm “ 2 ; the intrinsic, rest-frame, 2-10 keV 
luminosity is ~ 1.9 x 10 45 erg s -1 . Based on the optical and X-ray 
spectral properties (optical line widths, intrinsic N H and intrinsic 
2-10 keV luminosity in excess to 10 44 erg s -1 ) we classify this 
source as a Type2 QSO. 2XMMJ022256.9-024258 is detected in 
all the WISE bands (see Table 3). 

2XMMJ022256.9-024258 also has strong radio emission as¬ 
sociated with the radio source 4C-02.11. The total 1.4 GHz ra¬ 
dio flux from the FIRST survey l lWhite et afll 1 9971) is ~ 340 mJy 
corresponding to a radio power of ~ 1.7 x 10 34 erg s -1 Hz -1 at 
z= 1.004. The source is also detected in the lower spatial reso¬ 
lution NVSS survey dCondon et al.ll 19981) . The total NVSS radio 
flux density (power) is 429 ± 14 mJy (2.18 ± 0.07 x 10 34 erg s -1 
Hz -1 ); the radio flux density from the NVSS is a more reliable 
measure of the total radio flux, since the VLA-B configuration 
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used for the FIRST survey could miss some of the diffuse ex¬ 
tended emission. For this source we also found in the lit erature 
fluxes at 178 MHz (fmuH: = 2.6 ± 0.65 Jy jGowemnal 1 1 9671) 


and 4.85 GHz (/ 4 . 85GHz = 0.116 ± 0.012 Jv: iGriffith et~ 


19950 ; 


assuming a power-law model the derived radio spectral index is 
~ 0.94, consistent with the typical values of lobe-dominated AGN 
( iKellermann & Owenlll988ll . 

In Figure 3 we show the radio intensity contours from the FIRST 
radio survey overlaid on the optical image; the radio counterpart 
of 2XMMJ022256.9-024258 is a double radio source having the 
typical Fanaroff-Riley Type II morphology (i.e. sharp edge lobes 
and bright hot spots); the absolute optical magnitude ( M R - -24, 
assumi ng a K-correction of ~ 2 mag, typical of a late type galaxy, 
see lFukugita et al.ll 19951) and the total radio power (~ 2.2 x 10 34 
erg s' 1 Hz -1 ) are consistent with the FRII classification according 
to the dividing lin e in the M R — L,- n j lo plane b etween FRI and FRII 
radio galaxies fsee lGhisellini & CelottilEzOoH) . 

The projected separation of the two bright hot spots is — 41 arc- 
sec on the sky, corresponding to a physical projected size of ~ 0.33 
Mpc at the redshift of the source. Given the discussed radio prop¬ 
erties, the object is clearly a radio-loud and lobe-dominated AGN. 
In this respect, we note that the measured size and total radio lumi¬ 
nosity are fully consistent with that observe d in other radio quasars 
(see Figure 4 in lKuzmicz & Jamroz\ll2012i ). 

Considering the broad band properties discussed above, 
2XMMJ022256.9-024258 can be considered another example 
of the rare class of Radio Loud Type 2 QSO, with overall 
properties very s i milar to those shown by e.g. AXJ 0 843+2 942 
dDella Ceca etail 120031) 6C0 905+39 ferlund et al.1 120081) or 
4C+39.2 9 dGandhi et al.ll 20061). Indeed according to the data re¬ 
ported in lErlund et al. 00081 see their Figure 3) this object is one 
of the most powerful sources currently know; its intrinsic X-ray 
luminosity (~ 2 x 10 45 erg s' 1 ), if compared with its 178 MHz lu¬ 
minosity (~ 1.3 x 10 35 erg s' 1 Hz' 1 ), is more than a factor 10 above 
the locus sampled by the Narrow Line Radio Galaxies at z < 1.0 in 
the 3CRR catalog and similar to tha t observed in Bro ad L ine Radio 
Galaxies and QSOs (see Figure 3 in lHardcastle et al. f 20091) . 


• 2XMMJ 100038.9+050955; f x /f 0 =138 

The optical magnitude of the single object in the error circle 
(see Figure 2) is R=23.6. The gathered optical spectrum, unfor¬ 
tunately characterised by a very low SNR, seems to be rather 
flat and featureless, thus we have no redshift information for 
2XMMJ100038.9+050955. The source is detected in all the WISE 
bands while there is no radio detection at the source position. As¬ 
suming that 2XMMJ100038.9+050955 is an absorbed QSO and 
using t he relatio n between f x /f a and the intrinsic 2-10 keV lumi¬ 
nosity ferusa et al .20Kl see Section 5) we can estimate a redshift 
z~1.0 (in the range 0.6-1.6 taking into account a scatter of 0.5 dex 
in the relation, see Figure 4, left panel). At z=T.0 the intrinsic 2-10 
keV luminosity would be ~ 9 x 10 44 erg s' 1 while the intrinsic 
N h ~ 5 x 10 22 cm' 2 . 

• 2XMMJ 121026.5+392908; f x !f 0 =50 

This source was the only one not observe d during our GTC run; it 
is a wel l known BL Lac object at z=0.617 JCaccianiga et al .1120021 : 
IPlotkin et alj|20ldl . The observed optical magnitude (R=19.16), 
the radio (19 ± 0.7 mjy at 1.4 GHz from the NVSS) and the X- 
ray (~ 3.6 x 10' 12 erg cm' 2 s' 1 in the 2-10 keV band) fluxes 
allow us to compute multi-wavelength spectral indices El ( a RX ~ 
0.54; O'ox=0.61; or^o=0.50) which are fully consistent with the 


BL Lac classification dStocke et alJll99ll : ICaccianiga et alJlT999h . 
The measured a RX (< 0.8) suggests a high-frequency-peaked BL 
Lac object (HBL) in which the low-energy component of their 
usual, double-peak, SEP pe aks between the UV band and X-rays 
iPadovani & Giommilll995li . The observed X-ray spectrum is well 
described by a featureless power-law model (T = 2.23 4 ° q|) with in¬ 
trinsic N h = 6 x 10 20 cm ' 2 and intrinsic 2-10 keV luminosity of ~ 
6.4 x 10 45 erg s' 1 . In the WISE survey 2XMMJ121026.5+392908 
is detected only at 3.4 and 4.6 /rm (WISE band W1 and W2, re¬ 
spectively). 

. 2XMMJ121134.2+390054; f x /f 0 =50 

The optical object clearly visible at the centre of the 
error circle (see Figure 2) has a magnitude R=20.77. 
2XMMJ121134.2+390054 is a well known object since the 
epoch of the Einstein Extended Medium Sensitivity Survey 
dStocke et alj ll99lli and it is cla ssified as a BL Lac object 
(MS1209+3917: IRector et afll200d) . The observed optical mag¬ 
nitude, radio (10.6 ± 0.6 mjy at 1.4 GHz from NVSS) and X-ray 
(~ 8.2 x 10' 13 erg cm' 2 s' 1 ) fluxes imply multi-wavelength 
spectral indices (a RX = 0.58; a n X = 0.62; Q'go=0.56) fully consistent 
with t he BL Lac classification dStocke et al.ll99ll : ICaccianiga et al .1 
1 19991) . As for 2XMMJ121026.5+392908 discussed above, also 
2XMMJ121134.2+390054 can be classified as a HBL object. 

From the analysis of the optical spectrum shown in Figure 2 we 
reach the conclusion that the only clear feature is that at ~ 7048 A, 
most likely associated with [OIIj/13727 A at z=0.89; we estimate 
an observed equivalent width (EW) of ~ 12 - 15 A, that rescaled to 
z=0 imply an EW ~ 6 - 8 A, very close to the li mit used to classify 
a source as a BL Lac object (EW<5 A. IStocke et afll99lh . 

While our spectroscospic classification as a BL Lac object 
agrees with previous results, the redshift pro posed here (z=0.89 ) 
is significantly different from that proposed by IRector et al.1 12000 ) 
(z=0.602). These latter authors did not find any evidence of emis¬ 
sion line(s) in their spectrum, and based their tentative redshift de¬ 
termination on low SNR absorption features. We are confident that 
the line observed at ~ 7048 A is real and this line does not have any 
reliable identification if the object is at z=0.602. Moreover at the 
proposed redshift (z=0.89) we were able to reproduce quite well 
also the shape of the underlying optical continuum (see Figure 2). 
Finally we note that this source is strongly variable: at least a fac¬ 
tor of 5.5 in the X-ray domain (by comparing the measured 0.1-2.4 
keV flux with the flux in the ROSAT All Sky Survey) and a factor 
of 3 in the optical do main (by comparin g our magnitude with the 
V,„ as =20 reported in IRector et alj[2000[) . so this variability could 
e xplain the clear diff erent shape from the optical spectra reported 
in lRector et al.1 ( boOQ ) and that reported here (see Figure 2). 

At z=0.89 the observed X-ray spectrum is well described by a 
featureless power-law model (T = 2.2Uj?[?Z) with intrinsic N H = 
1.9 x 10 21 cm “ 2 ; the intrinsic 2-10 keV luminosity is ~ 3.7 x 10 45 
erg s' 1 . The intrinsic absorption is significantly in excess to the 
Galactic value and this is at odds with what is usual ly observ e d in 
BL Lacs and, in particular, in HBL objects (see e.g. lMassaro et al .1 
1201 ll and references therein). On the other end the X-ray spectrum 
could be intrinsically curved due to the bump expected in HBL, 
peaking between the UV and X-ray bands. We therefore fitted the 
X-ray data with a broken power-law model filtered by the Galactic 
absorption; we find a good fit (xl = 1.1 for 194 d.o.f.) with best¬ 
fitting parameters Ti = 1.57 4 q 4 , T 2 = 2.30 4 jj'J® and E c = 1.23 4 g 2 8 , 
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which are consistent with the results reported in iMassaro et all 
<201 lh on a sample of HBL sources. The intrinsic 2-10 keV lumi¬ 
nosity obtained with the broken power-law model is equal to that 
previously reported in Table 1. In analogy with the other confirmed 
BL Lac object in the sample (2XMMJ121026.5+392908), in the 
WISE survey this source is detected only at 3.4 and 4.6 fj m . 

• 2XMMJ123204.9+215254; /,//„ =1118 

Two optical sources (labeled with A and B in the finding chart 
reported in Figure 2) are clearly present inside the X-ray error cir¬ 
cle; they have an R magnitude of R=23.94 ± 0.41 (source A) and 
R=24.73 ± 0.64 (source B). An optical spectrum of source A, ob¬ 
tained summing up two independent exposures, is shown in Figure 
2 ; the spectrum has a clear red continuum and a significant line 
at ~ 6569 A. Assuming this line to be [OIIJT3727 A at z=0.763 
another observed feature is in very good agreement with being 
[OIII], 15007 A ; these two features are clearly seen in both the sin¬ 
gle exposures and are not associated with strong sky lines. We have 
also tried to take an optical spectrum of the source B, but the source 
is very faint and no useful information could be extracted from the 
very noisy spectrum. 

2XMMJ123204.9+215254 is detected in all the WISE bands dis¬ 
cussed here. The WISE position is coi ncident w ith a sou rce de¬ 
tected in the Kband (K= 18.07±0.05, see lDel Moro et al.l2009l l and 
both are significantly closer to the brightest optical source (source 
A), strongly suggesting that the most probable optical counter¬ 
part of the X-ray source 2XMMJ123204.9+215254 is the source 
A. In the following of this paper we therefore assume that the 
2XMMJ123204.9+215254 is spectroscopically identified with a 
Type 2 QSO at z=0.763. 

2XMMJ123 204.9+215254 was previously discussed in 
iDel Moro et alj <2009l) who report an infrared spectrum taken 
at the Subaru telescope (with the MOIRCS instrument) clearly 
revealing the presence of a line at the observed frame of 1.8837yum. 
Lacking optical spectroscopy and a deep optical image, these au¬ 
thors discussed several possible identifications for the infrared line 
and suggested the Ho- line at z=1.87 as its most probable origin. 
As discussed above our optical spectroscopy at the GTC suggests a 
lower redshift of z=0 .763. In this case the NIR feature reported in 
Pel Mor o et al.1 <20091 observed at 1.8837 yum) could be associated 
with the line complex Held 10 830+ Pay10941, a str ong and quite 
common feature in AGN (see iGlikman et al] l200dh . expected to 
be at ~ 1.92 /rm (at z=0.763), slightly higher than the observed 
feature. However the observed NIR line falls in a wavelength range 
dominated by the effect of atmospheric absorption, so it is difficult 
to determine the intrinsic line centroid; all in all, we consider the 
association of the infrared line with the HeLll0830+Payl0941 
complex as highly plausible. On the contrary none of the faint 
features that we see in the optical spectrum could be associated 
with any relevant emission lines from an AGN in the case of 
z=1.87. 

At the proposed lower redshift (z=0.763) the X-ray spectrum is 
described by an absorbed power-law model having a very flat pho¬ 
ton index, T = 1.31 + 0.23, and an intrinsic N H = 3. 6^2 x 10 22 
cm . The intrinsic 2-10 keV luminosity is ~ 1.9 x 10 45 erg s -1 
. Since the best-fitting photon index is rather flat we test the 
stability of the measured N H and luminosity assuming a typi¬ 
cal AGN photon index (T = 1.9); we derive an intrinsic N H = 
5-9-065 x 10 22 cm -2 and an intrinsic 2-10 keV luminosity of 
~ 2.3 x 10 45 erg s _1 in good agreement with the previous values. 
We note that 2XMMJ123204.9+215254 is not only the object with 
the highest f x !fo in this sample, but (to our knowledge) it is the 


non-transient spectroscopically identified source with the highest 
fxlfo discovered so far. 

• 2XMMJ135055.7+642857; f x /f„ =458 

A very faint object (R=25.0) is visible at the centre of the er¬ 
ror circle. Unfortunately the optical spectrum is very noisy and no 
features are clearly detected; we have no redshift information for 
this object. The X-ray spectrum is described by a power-law model 
with photon index T ~ 2 and absorbing column density greater 
than ~ 2 x 10 21 cm ~ 2 (obtained assuming z=0). An emission line 
is possibly detected in the EPIC-MOS spectrum (see the residuals 
in Figure 2, where the model does not include the line). However 
no compelling evidence of the presence of this line in the EPIC-pn 
is found, although different spectral binnings were tried. 

If this line is real, and associated with the Fe Ko- emission line 
(the most prominent feature in the X-ray spectrum of an AGN), the 
implied redshift would be z~ 0.65, the absorbing column density 
~ 4x 10 21 cm “ 2 and the intrinsic 2-10 keV luminosity ~ 2.5 x 10 44 
erg s _1 ; the line has an observed equivalent width of EW ~ 400 
eV. Interestingly the source is also a strong (and compact ) radio 
source, detected at 15 GHz (236 + 1 mJv. lRichards et al.l201lh . 8.4 
GHz (~ 376 mJv. lHealev etalll2007t) and at 1.4 GHz (183.5 + 5.5 
mJy, NVSS survey). The broad band spectral indices ( a RX = 0.88; 
a ox= 0.31; a^ 0 =1.15), computed with the fluxes reported above, 
and the possible presence of a Fe K» line in the X-ray spectrum 
are not consistent with a BL Lac classification. We remind that the 
X-ray source is point-like and detected almost in the center of the 
EPIC fields, so a high z cluster of galaxies is very unlikely. 

The radio spectrum is very flat (indeed almost inverted, with 
a maximum at ~ 8 GHz) suggesting that this source could be 
classified as a Giga-Hertz Peaked Spectrum radio source, a class 
of sources supposed to be young radio galaxies whic h are oft en 
characterised by high intrinsic absorption (see e.g. lGuainazzi et all 

hooch ). 

Assuming that this source is an absorbed QSO and using the 
relation between f x /f„ and the intrinsic 2-10 keV luminosity (see 
Section 5) we can estimate a redshift z~ 1.7 (z in the range 1.1 -2.7, 
inconsistent with the redshift estimated by the putative iron line), 
with an intrinsic 2-10 keV luminosity ~ 3 x 10 45 erg s _1 and an 
intrinsic Nh ~ 1.2 x 10 22 cm 2 . This source is detected in all the 
WISE bands. 

• 2XMMJ143623.8+631726; f x /f a = 55 

A faint object with optical magnitude R=22.16 is present in the 
centre of the error circle (see Figure 2). The optical spectrum of this 
object is shown in Figure 2; the two marked lines can be associated 
with [OII]T3728 and Hy+Hel (observed FWHM < 1100 kms^ 1 
) at z=0.893. The observed X-ray spectrum is described by an ab¬ 
sorbed power-law model with intrinsic N H = 1.46 x 10 22 cm ~ 2 ; the 
intrinsic 2-10 keV luminosity is ~ 8.4 x 10 44 erg s _1 . The source 
is detected in all the WISE bands; there are no radio detections at 
the source position. We classify this source as a Type 2 QSO. 


5 DISCUSSION 

Having defined a small but representative sample of bright EXO50 
objects, it is now instructive to compare their broad band proper¬ 
ties with those of other samples of AGN, both absorbed and unab¬ 
sorbed, from the literature. 
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5.1 The comparison sample 


We have assembled an heterogeneous sample of spectroscopically 
confirmed X-ray selected Type 1 and Type 2 AGN from a few sur¬ 
veys carried out during the last fe w years usin g XMM-Newton 
data: the COSMOS sur vey iB rusaetafl 2010h . the XMS sur- 
vev dBarc ons et all [20071 CarreraetjdJ, 2007) and the XBS survey 
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C accia niga ef aP l2008h . This first large 


comparison sample, composed by 882 Type 1 AGN and 487 Type 
2 AGN, will be used below (see Figure 4, left panel) to discuss 
the position of absorbed and unabsorbed AGN in the f x /f a vs. L x 
plane. 

As a second step, and in order to compare the optical-infrared 
photometric properties of the EXO50 objects with objects having 
similar luminosities and redshifts, we have considered in the com¬ 
parison sample only the high luminosity (intrinsic L ( 2 -io keV) > 10 44 
erg s _1 ) sources (i.e. the Type 1 QSO and Type 2 QSO) in the red- 
shift range between 0.6 and 1.1; we have cross-correlated this latter 
sample of sources with the WISE All-Sky source catalogue, using 
a searching radius equal to 2 arcsec (consistently with the results 
obtained between the cross-correlation of the EXO50 sources and 
the WISE catalog, see Section 3.3). 

The WISE-QSO sample to be used in the comparisons de¬ 
pends also on the WISE band(s) considered; to obtain meaning- 
ful colo ur distr ibutions we use a similar approach followed by 
iYan et all d2013l) and base our comparisons (see below) on the fol¬ 
lowing two WISE-X-ray selected QSO samples: the QSO Wl/2 
sample (S/N wl S* 7 and S /Nwi > 5; 94 Type 1 QSO and 11 Type 
2 QSO), the QSO W1/2/3 sample (S/N W i > 7, S/N W2 > 5 and 
S/N m > 3; 72 Type 1 QSO and 7 Type 2 QSO); the QSO W1/2/3 
sample is obviously a subset of the QSO Wl/2 sample 


5.2 The X-ray to optical flux ratio vs Lx and the optical-IR 
colour 

In Figure 4, left panel, we report the EXO50 objects discussed 
in this paper and the comparison sample of Type 1 and Type 

2 AGN from the literature. In the comparison sample 12 ob¬ 
jects have fxlfo > 50 and all but two □ have been classified 
as Type 2 QSO. Unfortunately the large majority of these 12 
sources (some of them will be highlighted in some figures re¬ 
ported below) are faint both in X-ray and in the optical domain, 
so only a couple of these interesting objects have been discussed 
in some detail in the literature. In particular we have the Type 2 
QSO XBS J021642.3-04 3553 (z=1.98; fJfo~ 200; L (Z - meV) ~ 

3 x 10 45 erg s _1 , see ISevergnini et all 120061 ) and the Type 2 
QSO XMMC_150.54703+1.61869 ( XID2028; z=1.59 ; f x /f 0 ~ 60; 
Lri ~ 8 x 10 44 erg s 1 , see lBrusa et al. I [20 1(1 iBrusa et al.l 
120151) ; both objects are considered as prototypes of the obscured 
QSO population at high redshift, where the coexistence between 
massive galaxies (~ few times 10 11 M & as estimated in the above 
papers) and powerful QSOs has been proved. Further more , in the 
case of XMMC_150.54703+1.61869, a recent paper dPerna et al .1 
l2014h reports the discovery of a massive outflow (> 630 M 0 yr 1 ), 
extending out to 10 Kpc from the central black hole and having a 
total energetic in full agreement with the prediction by AGN feed¬ 
back models. 

The dashed diagonal line in Figure 4 (left panel) corresponds 


6 These two sources, classified as broad line AGN, are CL0016.3 (z= 1.09) 
from the XMS survey and XMMC.150.13303+2.30324 (z=1.6) from the 
COSMOS survey. 


to the relation between f x lf 0 and the intrins ic 2-10 keV luminosity 
for obscured AGN initially proposed in lFiore et alj d 19991) and re¬ 


dly ] 

vised bv IBrusa et alJ |2010j), while the two dotted lines corresponds 
to a scatter of about 0.5 dex around this relation. As discussed in 
the previous section, 3 EXO50 objects are spectroscopically iden¬ 
tified as Type 2 QSOs (the sources number 1, 5 and 7 in the figure); 
taking into consideration the scatter that we see for the other Type 2 
AGN these 3 EXO50 objects seems to follow the same relationship 
between f x /f 0 and the intrinsic 2-10 keV luminosity. 

It is worth noting the position of the two EXO50 objects spec¬ 
troscopically identified as BL Lac objects (the sources number 3 
and 4); both objects seem to have a larger X-ray luminosity for a 
given f x lf 0 if compared with the absorbed AGN population. This 
is indeed expected, since BL Lacs have SED strongly dominated 
by the beamed emission from the jets. Their high f x /f a therefore 
is not due to obscuration (as for the Type 2 QSO where the optical 
luminosity is very likely dominated by the host galaxy) but to the 
intrinsic shape of the SED (see also below). 

In Figure 4 (right panel) we plot the f x jf D ratio as a function 
of the optical-infrared color (R-Wl) for our EXO50 objects and 
for the QSO W1/W2 sample defined in Section 5.1.The colour R- 
W1 can be considered a good proxy for the usual R-K colour and 
the vertical dashed line (R-W1=6.5) corresponds 0 roughly to R- 
K=5, the value used to define an object as a Extreme ly Red Object 
(ERO, see e.g. lElston et ali i 19881 iDaddi et al.ll200oh . X-ray emit¬ 
ting EROs studied so far, strongly suggest that the bulkof this popu- 
lation is composed by obscured AGN (se e e.g. lMignoli et al.l2004l 
ISevergnini et al.ll2005l Brusa et aPl2005l) and the results obtained 
here support this view[J 

A clear trend between the measured f x lf 0 and the optical- 
infrared colours is present. This trend, visible bot h for the Type 1 
QSO and for the Type 2 QSO samples (see also lBrusa et al.ll2005l) . 
could be partially explained as due to the effect of absorption (see 
also Section 5.6). As previously said, the nuclear optical fluxes are 
more depressed (indeed completely blocked in the case of Type 2 
QSO where we should see only the host galaxy) by circum-nuclear 
matter if compared with the infrar ed and X-ray fluxes. As already 
noted in ICaccianiga & Severgninil d201 lh the effect of absorption 
is important also in Type 1 AGN and this explain why we see the 
trend also in Type 1 QSOs. A further support to this hypothesis 
comes from the observation that, for the Type 1 QSOs in Figure 4 
(right panel) belonging to the XBS survey (for which we have X- 
ray and optical spectral information), there is a clear evidence that 
the sources with the higher f x /f 0 and redder colours are also those 
with the higher intrinsic Nu- 

Two out of three of the EXO50 sources classified as Type 2 
QSO are on the "EROs” side of the diagram (R-Wl > 6.5); the 
only object just outside the EROS locus is 2XMMI2256.9-024258 
(the EXO50 object at z~l associated with the double and bright 
radio source 4C-02.il), although its position in the plot is not sig¬ 
nificantly different from the other sources considering the scatter 
on the f x lf 0 vs R-Wl relation. The position of the two unidenti¬ 
fied EXO50 in the f x lf 0 vs. R-W1 plot (Figure 4, right panel, open 
squares), in the R-W2 vs. W1-W2 plot (Figure 5, open squares) 
and in the f x !fo vs. R-W2 plot (Figure 7, open squares) strongly 


7 This correspondence has been derived from the analysis (not reported 
here) of the sources in the COSMOS survey where we have both R, K and 
W1 magnitudes. 

8 We note that both objects from literature discussed in the previous sec¬ 
tion are EROs (R-K>5 for XBS J021642.3-043553 and R-K=6.46 for 
XMMC.150.54703+1.61869). 
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Figure 4. Left panel', relationship between the measured f x / f„ (in logarithmic units) and the intrinsic 2-10 keV luminosity for our sample of EXO50 objects 
and for the comparison sample of AGN assembled as described in Section 5.1. Large squares: our sample of EXO50 objects (filled-red for Type 2 QSO, 
open-red for unidentified objects and open-crossed-green for BL Lac objects) the numbers close to the large squares mark the EXO50 objects as reported in 
Table 1 and Table 3. The labels a) and b) mark the position of the objects XBSJ021642.3-043553 (z~ 2, source a) and XID2028 (z~ 1.6, source b) discussed in 
section 5.2. Filled red circles: Type 2 AGN from the comparison sample; Small blue crosses: Type 1 AGN from the comparison sample; The dashed diagonal 
corresponds to the relation between f x /f a and the intrinsic 2-10 keV luminosity for obscured AGN, while the two dotted lines corresponds to a scatter of about 
0.5 dex around this relation. Right panel', the f x lf a ratio as a function of the optical-infrared color R-Wl (a proxy of the usual R-K colour) for the EXO50 
objects and for the sample of QSO W1/W2 as described in the text; the vertical dashed line (R-W 1 =6.5) corresponds roughly to R-K=5. The horizontal dashed 
line correspond to fxlfo = 50. Symbols and colours are as in the left panel. 


suggest an obscured QSO nature; if we assume that the two uniden¬ 
tified sources follow the f x /f 0 —L x relation, their intrinsic L x would 
be 2.5x10 44 -2.5x10 45 erg s _1 (which would imply aredshift in the 
range 0.6-1.6) for 2XMMJ100038.9+050955 and9xl0 44 -9xl0 45 
erg s _1 (which would imply a redshift in the range 1.1 — 2.7) for 
2XMMJ135055.7+642857. 


5.3 The optical-IR flux ratio vs the IR colours 


The optical-infrared colours (R-W2) for our EXO50 objects and 
for the comparison QSO W1/2 sample as a function of the infrared 
colour W1-W2 are shown in Figure 5. We have also marked the 
place (Wl-W2>0.8 and R-W2> 6) where, according to 1 Yan et all 
d2013h , it is possible to select Type 2 QSO candidates. The two 
EXO50 objects spectroscopically identified as BL Lac objects seem 
to have optical-IR p roperties different fr om the bulk of the AGN 
population (see also iMassaro et al.ll20fll for similar results based 
on WISE MIR colours). We stress that our EXO50-Type 2 QSO 
objects are amongst the most extreme sources in these diagrams 
(see Figure 15 in Yan et al.ll2013l for a comparison with the SDSS 
QSO sample). 

As can be seen, all the QSOs, regardless of their spectro¬ 
scopic or f x /f„ properties, span a very similar range of the infrared 
(W1-W2) colour, although there are a couple of our EXO50 ob¬ 
jects with extreme W1-W2 colours (source 2 and source 5, namely 
2XMMJ 100038.9+050955 and 2XMMJ123204.9+215254). On 
the contrary the optical-infrared colours (R-W2 or R-Wl, see also 
Figure 4, right panel) have to be in some way correlated with the 
X-ray to optical flux ratio properties, and with the optical spectro¬ 
scopic classification (see also the results reported in Figure 7). 


The observed quite similar range of the W1-W2 colours in 
Type 1 and Type 2 Q SOs (i n agreement with similar results previ¬ 
ously reported in e.g. lMateos et al.lF2013h suggests that the absorp¬ 
tion is probably not the only (and the more important) parameter in 
shaping the AGN infrared continuum (between ~1.7 and 3/rm rest 
frame assuming an average redshift of 0.8). On the contrary, the 
observed separation of the different classes of AGN in the R-W2 
colour clearly indicates that the absorption could be an important 
parameter in shaping the optical-infared observed properties. We 
will discuss these topics in Section 5.6. 


5.4 The IR colours diagnostic plot 

In Figure 6 we show the WISE MIR colours (this time as a ratio 
between fluxes) for our EXO50 objects and for the QSO W1/2/3 
sample defined above; the two EXO50 BL Lac objects are not re¬ 
ported here, since they are not detected at 12/um. The solid wedge 
represents the locu s populated by luminous AGN ac cording to 
IMateos et alJ d2012l) . As shown in iMateos et al.l d2012l) and con¬ 
firmed here, we see that the technique is highly effective at identi¬ 
fying both Type 1 and Type 2 QSOs. The position of the different 
kinds of QSOs in the WISE MIR colour plot does not seem to sug¬ 
gest a big difference between absorbed and unabsorbed AGN in 
the MIR domain, in good agreement with recent results reported in 
IMateos et alJd2013l) . 

It is worth noting that three EXO50 objects clearly stand 
out from the zone populated by the bulk of AGN; two of 
these objects are spectroscopically identified and, interestingly, 
are the sources with the highest X-ray luminosity (namely 
2XMMJ123204.9+215254 and 2XMMJ022256.9-024258). 
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Figure 5. The optical-infrared R-W2 colour vs. the infrared colour Wl- 
W2 for the EXO50 objects and for the comparison QSO Wl/2 sample (see 
Section 5). The zone in the plot enclosed by the dashed line (R-W2>6 and 
Wl-W2>0.8) co rresponds to the place where it is possible to select Type2 
QSO candidates lYan et al]|2013l) . Symbols and colours are as in Figure 4; 
we have also marked the sources with f x lfo > 50 in the comparison QSO 
samples using empty squares. The numbers mark our EXO50 objects as 
reported in Table 1 and Table 3. The labels a) and b) mark the position of 
the objects XBSJ021642.3-043553 (z~ 2, source a) and XID2028 (z~ 1.6, 
source b) discussed in section 5.2. 


Figure 6. WISE MIR colours for the EXO50 objects and for the compari¬ 
son QSO Wl/2/3 sample (a subset of the QSO W1/W2 sample reported in 
Figure 5, see Section 5.1 for details). Note that the two BL Lacs are not 
reported here since they are not detected at 12 /im. The solid wedge repre¬ 
sents the AGN locus according to lMateos et alj20~12l (see Section 5 details). 
Symbols and colours are as in Figure 4; we have also marked the sources 
with fxl fo > 50 in the comparison QSO samples using empty squares. The 
numbers mark the EXO50 objects as reported in Table 1 and Table 3. The 
labels a) and b) mark the position of the objects XBSJ021642.3-043553 
(z~ 2, source a) and XID2028 (z~ 1.6, source b) discussed in section 5.2. 


5.5 The relation between EXO50 objects and EDOGs 

The three EXO50 sources on the top-right zone of Figure 5 (the 
sources 2, 5 and 6) have R-W4 above 15 (see Table 1 and 3), 
implying a / 24 ^ m ///? well in excess of ~ 3000. These 3 sources 
can be classified as Extreme Dust Ob scured Galaxies (EDOGs; 
fiAumlfR > 2000; lLanzuisi et al]|2009h , i.e. extreme ex amples of 
Dust Obscured Galaxies (DOG, fa > 1000 see I Dev et ail 
l2008LiFiore etai]|2008h . a class of galaxies firstly discovered with 
the Spitzer Space Telescope. In these objects, the high infrared-to- 
optical ratios imply that large amounts of dust are absorbing the op¬ 
tical light and re-emitting it in the infrared. DOGs and EDOGs may 
play an important role in the formation of the most massive galaxies 
in the local Universe and, in particular, represent an i mportant evo- 
lutionary step in the AGN-galaxy connection (see lLanzuisi et al.l 
120091 and reference therein). 

A sample of 44 EDOGs from the SWIRE survey having X-ray 
coverage (with XMM-Ne wton or Chandra) has been selected and 
studied bv lLanzuisi et alJd2009f) ; the source redshifts (both spectro¬ 
scopic and photometric), are in the range between ~0.7 and ~2.5. 
About 95% of the detected sources 0 (23 objects) are consistent 
with being obscured by neutral gas with an intrinsic column density 
(mostly derived using X-ray hardness ratios) typically in the range 
between 10 22 to few times 10 23 cm' 2 . Their intrinsic 2-10 keV lu¬ 
minosities (in excess to 10 43 erg s' 1 ) fall well within the AGN X- 
ray luminosity range and ~ 55% of them can be classified as Type 2 

9 An inadequate X-ray coverage is probably the cause of the non-detection 
of the remaining ones. 


quasars, on the basis of their absorption properties and X-ray lumi¬ 
nosity. It is also worth noting that, among the EXO50 sources de¬ 
fined here, about 40% can be classified as EDOGs (~ 60% exclud¬ 
ing the BL Lac objects). Unfortunately, two of our EXO50 EDOGs 
are still unidenti fied, prevent ing us to make a proper comparison 
with the work of lLanzuisi et all J20091 ). We can only say that the 
absorption properties of the only identified EXO50 EDOG (source 
number 5) are fully consistent with the main results obtained from 
lLanzuisi et al.1 d2009h : the source is an absorbed QSO with an in¬ 
trinsic absorbing column density ~ 3x 10 22 cnT 2 . For the remaining 
two the Type 2 QSO hypothesis is very likely. 

5.6 Are the EXO50 Type 2 QSO objects different from 
standard X-ray selected Type 2 QSO? 

Overall the EXO50 Type 2 QSO do not seem to be different from 
standard X-ray selected Type 2 QSOs in terms of nuclear absorp¬ 
tion (in the range of few times 10 22 cm' 2 up to 8 x 10 23 cm' 2 ), so 
other factors may play a possible role in explaining their extreme 
properties. 

To investigate how the f x /f 0 and the optical-to-infrared 
colours change as a function of a) the intrinsic absorption of the 
AGN component; b) the luminosity contrast between the intrinsic 
AGN and the host galaxy; and c) the redshift, we have used a very 
simple toy model to describe the SED of the AGN+host galaxy 
complex. 

As unobscured AGN SED we have used the Type 1 QSO UV- 
to-IR template described in lPolletta et al.l i2007l) adding the X-ray 
component assuming a f x jf„ (as defined in Section 2) =s 1 (typi- 
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Figure 7. The f x lfo ratio as a function of the optical-infrared color for our 
EXO50 Type 2 QSO objects and for the comparison sample. We also report 
the expectations from the toy model at z=0.6 and z=l.l for two values of 
optical obscuration (Av=5.27 and 52.7 corresponding to Nh — 10 22 and 
10 23 cm -2 , respectively) and three values (10, 20, 100) of AGN/host ratio 
at ~ 4.6 fim ; the big star marks the position where the toy model converges 
for all the AGN/host ratio used here when Nh < few times 10 20 cm -2 (i.e. 
at the unabsorbed Type 1 view). Symbols and colours are as in Figure 4. 
The numbers mark our EXO50 objects as reported in Table 1 and Table 3. 
The labels a) and b) mark the position of the objects XBSJ021642.3-043553 
(z~ 2, source a) and XID2028 (z~ 1.6, source b) discussed in section 5.2. 
See Section 5.6 for details on the modelling. 


cally observed in unabsorbed AGN), while for the host galaxy we 
have used a Sc galaxy template; similar trends to those discussed 
below are obtained if we assume a different morphology for the host 
galaxy template. Full details on the AGN and host SED as well as 
on the procedure summarised below can be found in IBalio et all 
J2014h . 

The UV-to-IR AGN template is absorbed at th e source r edshift 
using the extinction curve of the Galactic center dChiar & Tielend 
1200611 before summing it up with the galaxy template; the X-ray 
absorption is tied to the optical obscuration assumin g a Galacti c 
gas to dust ratio (N H /Av = 1.9 x 10 21 cm -2 . iBohlin et al.lll978i l. 
The variables of this simple toy model are thus the redshift, the 
AGN dust extinction (and consequently the X-ray absorption) and 
the ratio between the intrinsic AGN emission and the host galaxy 
emission (normalised for our convenience at ~ 4.6 gm, observed 
frame). 

In Figure 7 we show the expectations of this toy model in the 
fxlfo vs. R-W2 plane at a redshift of 0.6 (solid line) and 1.1 (dashed 
line, respectively the lower and the upper envelope of the investi¬ 
gated redshift range) for two values of dust extinction (Av=5.27 and 
52.7, corresponding to N H — 10 22 and 10 23 cnT 2 , respectively) and 
for three values of the AGN to host galaxy ratio (10, 20 and 100). In 
the figure we have also reported the position of the sources belong¬ 
ing to the QSO Wl/2 sample as well as our Type 2 QSO EXO50 
sources. The big star marks the position where the toy model con¬ 
verges, for all the AGN/host ratios used here, when Nh < few times 
10 20 cnT 2 . As can be seen, the expected position of the unabsorbed 


AGN is fully consistent with the data for the Type 1 QSO in the 
comparison sample. 

Figure 7 clearly suggests that several factors may concur in 
explaining the extreme properties of the EXO50 Type 2 QSO. For 
a fixed absorption or AGN/host ratio, the higher is the redshift, 
the higher i s the f y /f„ and the redder are the optical-to-infrared 
colours (cf. iFiore et al.ll20031 ). However from a deeper look to the 
data, this redshift effect alone is not able to explain the extreme 
properties of all the EXO50 objects; for example the source num¬ 
ber 5 (namely 2XMMJ123204.9+215254 at z=0.76) can not be ex¬ 
plained with the redshift effect but requires very high AGN/host 
ratios; a very high AGN/host galaxy ratio indeed may be required 
also for 2XMMJ 100038.9+050955 and 2XMMJ135055.7+642857 
(sources 2 and 6, respectively) if their redshifts turn out to be sim¬ 
ilar to that of 2XMMJ123204.9+215254. For comparison we note 
that the 8 Type 2 QSO with f x /f D in the range between 5 and 15 
reported in Figure 7 belong to the XBS survey and have been stud¬ 
ied in detail bv lBallo et al J d2014l) ; these sources have an AGN to 
host galaxy ratio at 4.6 fim (observed frame) between a few and = 
20 . 

A wide range in the intrinsic luminosity ratio between the 
AGN and the host galaxy could also explain the observed wide 
range of the W1-W2 colours in Type 1 QSO and Type 2 QSO 
(see Section 5.3), with the AGN component increasing, on aver¬ 
age, going from the left (Wl-W2=0, corresponding to a ^ -2) to 
the right (Wl-W2=2, corresponding to a =4.3) of the abscissa in 
Figure 5. Finally, as already discussed in the previous section, the 
sources with the highest f x /f D are also the sources with extreme 
optical-IR properties. Therefore a high AGN/host ratio (along with 
the large amount of dust absorption) could be a natural explanation 
to account for the extreme properties for a part of the EDOG popu¬ 
lation. This also implies that the infrared luminosities measured in 
these objects are mainly due to the AGN component rather than to 
starburst activity. 

5.7 The BL Lacs in the EXO50 sample 

Finally, although the statistics of EXO50 objects is very limited, 
it is worth noting that ~ 30 per cent of the EXO50 sources in the 
bright flux regime explored by our survey is represented by BL Lac 
objects. At the 2-10 keV flux limit of ~ 1.5xl0 -13 erg cnT 2 s -1 the 
derived surface density is ~ 3 x 10 -2 deg -2 . This surface density is 
a lower limit to the real surface density of BL Lac objects, since we 
are considering here only the BL Lacs with f x / f 0 > 50, and corre¬ 
sponds to a fraction between ~ 8% and ~ 30% of the density of X- 
ray selected BL Lacs at similar fluxes (the exac t fraction dependi ng 
on their cosmological evolution properties, e.g. lWolter et al.ll99ll ). 
Their position in the several plots discussed here deviates signifi¬ 
cantly from the rest of the sample. The very high f x lf 0 selection 
used here allows us to pick up the most extreme BL Lacs (HBL 
type) in the a ox - a rx plane (see also lCostamante et alj200ll) . This 
is clearly evident in Figure 8, where we show the a n}( -a R x pl ane for 
the BL Lacs objects in the BZCAT catalog dMassaro~et aP l2009l) 
with the position of the two EXO50 sources identified as BL Lacs 
overlaid; only less than 1 percent of the BL Lacs discovered so far 
have broad band spectral indices as extreme as the ones discovered 
here using the f x l f 0 > 50 selection. 

BL Lac objects are sources emitting non-thermal radiation 
across the entire electromagnetic spectrum from a relativistic jet 
that is viewed closely along t he line of sight, thus ca using strong 
relativistic amplification (e.g. lUrrv & Padovanill 19951) . Their SED 
is characterised by two broad peaks, now almost commonly in- 
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Figure 8. The ctox-^Rx plane for the BL Lacs objects in the BZCAT catalog 
(black open circle) and for the two EXO50 sources (open-crosses-green 
symbols) spectroscopically identified as BL Lacs reported in this paper 

terpreted as due to synchrotron and inverse Compton radiation. 
In this respect, the BL Lac objects reported here should have 
the synchrotron peak e mission around 0.1 keV (see Figure 1 in 
iCostamante et alj 120011 ). indicating the presence of high energy 
electrons: these extreme BL Lacs are good candidates for TeV 
emission that, if detected at these redshifts, could provide important 
constraints both on the physical process at work in such sources, 
and on the intensity and spectru m of the diffuse Extragalactic Back¬ 
ground Light (see e.g. Icostamantell2013l) . Presently (end of July 
2014) there are about 50 blazars detected at TeV energies 0. Two 
out of the 6 detected blazars at z> 0.3 are HBL (with z=0.34 and 
z>0.6) like the ones discussed here; the X-ray and radio fluxes 
of these two high-z TeV detected HBL are comparable with the 
fluxes of our sources and, therefore, we expect that these objects 
are potentially detectable at TeV energies even with the facilities 
currently available. 


6 SUMMARY AND CONCLUSIONS 

Obscured QSOs are expected to display large values of X-ray to 
optical flux ratio, fxlfo , where f x refers to the observed 2-10 
keV flux (corrected for Galactic absorption) and f„ is the optical 
flux derived from the R magnitude. In this paper we have defined 
a small, but statistically complete and representative sample of 7 
point-like X-ray sources characterised by an extreme X-ray to opti¬ 
cal flux ratio, f x /f D >50 (i.e more than 15 times the average values 
of unobscured broad line and soft X-ray selected AGN), and with 
f x > 1.5 x 1CL 13 erg cnr 2 s _1 . The sky coverage investigated to 
find out these 7 objects is about 60.4 sq deg, implying a density of 
~ 1.2 EXO50 sources/deg 2 at the sampled fluxes. We have spec¬ 
troscopically identified (using the GTC data as well data from the 
literature) 5 EXO50 objects; two sources remain still unidentified, 

10 see http://www.asdc.asi.it/tgevcaf/index.php 


although the data at other wavelengths discussed here strongly sug¬ 
gest an obscured QSO nature. 

We have discussed in detail the broad band properties of each 
EXO50 source in our sample, and we have compared them with 
those of a comparison sample of X-ray selected AGN/QSO from 
the literature. 

The main results of this investigation are the following: 

a) about 70 percent of the EXO50 sources in the bright flux 
regime explored by our survey are associated with obscured AGN 
(N h > 10 22 cnr 2 ), spanning a redshift range between 0.75 and 1 
and characterised by 2-10 keV intrinsic luminosities in the QSO 
regime (e.g. well in excess to 10 44 erg s -1 ). This result confirms 
the suggestion that the f x /f a ratio can be used as a proxy of ob¬ 
scuring material and it is a ver y efficient tool for the selection of 
obscured QSOs (see also lLanzuisi et al1l2013l) . If compared with 
other samples of Type 1 and Type 2 QSO, the EXO50- Type 2 QSO 
objects are amongst the most extreme sources in the several optical- 
infrared diagrams investigated. Overall the EXO50 Type 2 QSO do 
not seem to be different from standard X-ray selected Type 2 QSOs 
in terms of nuclear absorption; a very high AGN/host galaxy ratio 
seems to play a major role in explaining their extreme properties. 
Interestingly, 3 out of 5 EXO50 Type 2 QSO objects can be clas¬ 
sified as EDOGs {fiAiiml /r > 2000), a source population that may 
play an important role in the formation of the most massive galax¬ 
ies in the local Universe. Similarly to what we have found for the 
EXO50 Type 2 QSOs, we suggest that a very high AGN/host ratios 
(along with the large amount of dust absorption) could be the nat¬ 
ural explanation for the obser ved properties for part of the EDOG 
population. Two recent papers (Brusa et al.l2015l .lPerna et al .l2014l) 
find the presence of massive outflows in objects similar to the ones 
discussed here. Unfortunately the optical spectra presented in this 
paper does not have the quality required to study in any detail the 
presence of outflows; an approved LBT programme on this topic is 
ongoing; 

b) the remaining EXO50 sources are represented by BL Lac 
objects. Interestingly the very high f x /f„ selection used here allows 
us to pick up rather extreme BL Lacs, which are good candidates 
for TeV emission. If detected, they could provide important con¬ 
straints both on the physical process at work, and on the intensity 
and spectrum of the diffuse Extragalactic Background Light. 
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